Motion of an impurity particle in an ultracold quasi-one-dimensional gas of hard-core 

bosons 
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The low-lying energy eigenstates of a one-dimensional (ID) system of many impenetrable point 
bosons and one moving impurity particle with repulsive zero-range impurity-boson interaction are 
found for all values of the impurity-boson mass ratio and coupling constant. The moving entity 
is a polaron-like composite object consisting of the impurity clothed by a co-moving gray soliton. 
The special case with impurity-boson interaction of point hard-core form and impurity-boson mass 
ratio nii/m unity is first solved exactly as a special case of a previous Fermi-Bose (FB) mapping 
treatment of soluble ID Bose- Fermi mixture problems. Then a more general treatment is given using 
second quantization for the bosons and the second-quantized form of the FB mapping, eliminating 
the impurity degrees of freedom by a Lee-Low-Pines canonical transformation. This yields the exact 
ground state (total linear momentum q — 0) and exact boson-impurity distribution function in the 
thermodynamic limit for arbitrary mi/m and arbitrary impurity-boson interaction strength. These 
results are then extended to states with g > 0. 

PACS numbers: 03.75.-b, 67.85.-d 
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Due to the rapidly increasing sophistication of exper- 
imental techniques for probing ultracold gases, theoret- 
ical emphasis has shifted from effective field approaches 
to more refined methods capable of dealing with corre- 
lations. When confined in a de Broglie waveguide with 
transverse trapping so tight that longitudinal energies 
are less than the transverse vibrational excitation energy 
quantum huj±, an ultracold Bose gas becomes effectively 
one-dimensional (ID) with a confinement-induced res- 
onance in the ID scattering length [l|. In units such 
that h = 1, the dimensionless ID coupling constant is 
jB — it^'Qid/''^ where gfj^ is the coupling constant for 
zero-range interactions gfj)S{xj — xi) dl of Lieb-Liniger 
(LL) form [2], and n is the ID density n = N / L with N 
the number of bosons and L the cell length for periodic 
boundary conditions. At low densities where 75 ^ 1, the 
boson-boson interaction reduces to impenetrable point 
form, the so-called Tonks- Girardeau (TG) limit for which 
the exact many-body ground and excited states were 
found some 48 years ago by the Fermi-Bose (FB) map- 
ping method [3]. The fermionization predicted therein 
has recently been experimentally confirmed \% |5| , lead- 
ing to much recent theoretical and experimental activity 
on various properties of such strongly-correlated quasi- 
ID ultracold gases. It was found recently [6] that impu- 
rity atoms in an ultracold Bose gas with velocity Vi have 
collision cross sections with the bosons which increase 
dramatically when Vi > c where c is the speed of sound 
in the Bose gas. Here we shall investigate the behavior 
of a moving impurity particle in an ultracold ID gas of 
impenetrable point bosons, assuming a repulsive, zero- 
range impurity-boson interaction of arbitrary strength, 
and arbitrary impurity-boson mass ratio. 



Hard-core impurity in a TG gas: Consider first the 
case Np = 1 of our previous exact solution of the prob- 
lem of N hard-core bosons and Np hard-core fermions in 
ID [7]. One starts from a "model wavefunction" '^m — 
J2p £{P)ui{Pxi) ■ • ■ UN{PxN)uN+i{y) where y is the im- 
purity position and the sum runs over all (A^+1)! possible 
permutations of these variables including permutations 
exchanging bosons with the impurity, s{P) is the usual 
±1 sign of the permutation, and ui, ■ • • , uat+i are A^ + 1 
orthonormal orbitals occupied by all A^ bosons plus the 
one impurity. $jv/ vanishes at all the points Xj — xg and 
Xj — y required by the hard-core constraints, and its im- 
proper symmetry under boson-boson and boson-impurity 
exchange can be repaired by a generalized FB mappin, 
function A = ni<j«?<jv sgn(a;j -x^) H^^i sgn(xj -y) 
where the sign function sgn(a::) is -|-1 (—1) if a; > (a; < 
0). Then the physical wavefunction is ^(xi, • • • , xm] y) = 
A{xi,- ■ ■ ,XN',y)'^ m{xi,- ■ ■ ,xi^;y). The model wave 
function ^1 m is an exact many-body energy eigenstate 
if the orbitals u,y{x) are eigenfuctions of the single- free- 
particle Hamiltonian —q^^t with eigenvalues e^, and 
the physical state ^ = A'^ m is an energy eigenstate with 
eigenvalue X^i/^i'- ^^ there is no external potential and 
the system is on a ring of circumference L with the Xj 
and y measured circumferentially, the u^ are plane waves 
j^-i/2^ik„x ^i^]^ 1^^ which are integer (-l-,-,0) multiples of 

^ and e^ = 7^. 

L ^ 2m 

Assuming A^ to be even, the model ground state ^om is 
a filled Fermi sea of the A^ -|- 1 lowest plane-wave orbitals 
with —kp < k^ < kp where kp = rnr and n = N/L, i.e.. 
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■^ and has total linear 



(or angular, on a ring) momentum zero. For this moving 
impurity problem, one needs excited states of nonzero to- 




FIG. 1: Excitation energy branches as functions of the mo- 
mentum transfer q for Ai' + 1 = 31 particles. The bold solid 
lines indicate the lowest-energy (hole) excitations, the thin- 
ner solid lines indicate the lowest excitations of both Fermi- 
surface (FS) and umklapp types, the dashed lines are the pro- 
longations of the FS and umklapp excitations to higher-energy 
values, and the underlying dashed parabola is q^ /2{N + l)m. 



tal momentuin q. The energy landscape is identical with 
that of A^-|- 1 TG bosons; only the physical interpretation 
is different. These excitations are the jb ^ +oo limit 
of the types I and II excitations of the LL gas, Fig. 4, 
p. 1620 of [3|, and are shown in Fig[TJ The lowest exci- 
tations, indicated by the bold solid lines, correspond to 
the type II excitations of [2'| and correspond to promot- 
ing a particle from kp — q + 2tx jL to hp + 2n / L or from 
—kp + \q\ — 2-k/ L to —kp — 27r/L, thus leaving a hole in 
the interior of the Fermi sea. As q increases from zero to 
2k p the hole moves to the left from kp to —kp, at which 
point the whole Fermi sea has been shifted to the right by 
27r/L. Here c = nh/2m is the speed of sound in the TG 
gas [3|. The excitation energy is periodic in the center 
of mass system with period 2mc = 2kp [9|, and the solid 
lines for |g| > 2mc correspond to this periodicity, the un- 
derlying dashed parabola q^ /2{N + l)m being the result 
of transformation from the center of mass to the labora- 
tory system, and the bold sohd hues correspond to hole 
excitations from the original and displaced Fermi surfaces 
kp + j2tt/L and -kp + j27r/L with j = 0, ±1, ±2, • • • . 
When j reaches {N + l)/2 the flow speed Vj — j2-K/mL 
reaches the speed of sound c, the slope in the entire inter- 
val qj-i < q < qj becomes > so that the energy barrier 
disappears, metastability is lost, and Mach wave drag 
commences. The boson-impurity distribution function 
Pbiix-y) = N{N-1) J\^{x,X2,- ■■ ,XN;y)\'^dx2 ■■ -dxN 
is found in the thermodynamic limit to be pbi{x — y) — 
n?[l—JQ{kf{x — y))\ where jo(C) = sin ^/^ is the spherical 
Bessel function of order zero. This is independent of the 
flow speed Vj and agrees with [Sj] . 

Quantized field representation: Generalize now to the 
case where impurity and boson masses are unequal and 
boson- impurity interactions are of the form \iji5{xj — 



y). The impenetrable point boson-boson interaction 
will eventually be used, but initially assume it to be 
\i,5{xj — x^). The method used is motivated by a pre- 
vious 3D weak-coupling treatment [8|. We use second 
quantization for the bosons while retaining Schrodinger 
representation for the impurity, with position variable y 
and momentum operator pi — --§-■ Then the Hamilto- 
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where tpB{x) and ipg{x) are the Bose field annihilation 

and creation operators, psiy) = V'b(2/)''/'s(2/) is the boson 
density operator at the impurity position, and we assume 
units wherein the boson mass m = 1. Now enlarge the 
state space to also contain Fock states with Fermi field 
creation operators 4'p{x) where 4'f and 4'p satisfy the 
usual anticommutation relations, but commute with the 
Bose fields, and both Bose and Fermi fields commute with 
the impurity position and momentum operators. On this 
enlarged Hilbcrt space, introduce a canonical transfor- 
mation which is the second-quantized form [10,] of the 
FB mapping [3|: 
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Then Hb is transformed into 
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Hp = I dx i/jp 
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where pp{y) — i^p{y)4'F{y) is the fermion density op- 
erator at the impurity position. The zero-range boson- 
boson interaction term in Hb has disappeared because 
its transform involves i^p{x) and its hermitian conjugate, 
which vanish by the anticommutation relations. Because 
of nonlocality of the FB mapping transformation, the 
transform of the boson kinetic energy operator does not 
have the simple form above on the whole Hilbert space, 
but on the subspace of states \^b) satisfying the im- 
penetrable point hard core constraint ^/'^(x)|<I'b) = it 
reduces to this form [3J, llO| . 

Next we make a second canonical transformation of the 
type used by Lee, Low, and Pines (LLP) in their theory 
of the polaron ll| which takes advantage of total linear 
momentum conservation so_as to eliminate the impurity 
dynamical variables, as in ^S]. It is effected by a unitary 
operator U = e"*^^^ where pp = J dx 'il)p{x)j-^ipp{x) 
is the fermion momentum operator, with the result 
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Then ^ is transformed into 
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where Tp, the fermion kinetic energy operator, is the first 
term of Eq. ([3]). Since y is absent from 7i, pi commutes 
with Ti. and may be replaced by its eigenvalue q. Any 
energy eigenstate of Ti. with eigenvalue E can be written 
as a direct product |$) = L~^/^e'*^|$F) where |$f) is 
independent of y. Then defining |$^) — f7|$) one finds 
that Hpl^g) = Egl^g) and P\^q) = q\^q) provided that 
Tlql^q) = Eql^q), whcrc Eq IS the energy eigenvalue of 
the given direct product eigenstate of Ti, P = Pi + pf 
is the conserved total linear momentum (impurity plus 
fermions), and Tiq is the Hamiltonian obtained from Ti. 
by replacing the impurity momentum operator pi by the 
c-number q. With this replacement Eq. ([5]) becomes 
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After mapping back to the physical Hilbert space (bosons 
plus impurity), q is found to be the total linear momen- 
tum in the laboratory frame. The ground state has g = 0, 
in which case two of the last three terms in © vanish 
identically. We will consider that case now and will gen- 
eralize later to the case q > 0. For q = Q the mean 
fermion momentum {pp) vanishes by symmetry, so the 
expectation value of the last term in ^ is proportional 
to the mean square fluctuation oi pp. This fluctuation 
term will be shown later to be negligible in the thermo- 
dynamic limit, and will be dropped now. 

Orbitals: With the last two terms dropped and the zero 
of energy shifted by q'^/2mi, Eq. © is the Hamiltonian of 
free fermions in the field of a potential Xi,i5{x) centered 
on the impurity. Its iV-particle ground state is a filled 
Fermi sea of the lowest N orbitals 4)k{x) of the single- 
particle Schrodinger equation H(j)k{x) — ek4'k{x) with 
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where ^pd- is an L-periodicized delta function, 6pcr{x) = 
X]^-oo S{x + i^L). Our approach here is a generalization 
of a previous treatment of dark solitons in a TG gas [I4I . 
As in ^j we assume that N is odd so that the mapping 
back to the boson Hilbert space will be periodic. It is 
convenient to choose real eigenstates of ([7]) since the A^- 
fermion ground state has zero total linear momentum, 
and this will be satisfied automatically if all the orbitals 
have zero flux. We generalize Eq. (6) of 12] to 



'Pk^\^) = A4-[sin(fc|x|) + AfcCOs(fca:)] 
'(a;) = ^/2/Lsm{kx) 
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where the allowed /c-values for both the even-parity or- 
bitals (pf, and the odd-parity orbitals </)[. are k — 
2tt/L, Att/L, ■ ■ ■ . These orbitals satisfy periodic bound- 
ary conditions on the interval —L < x < L, and are to be 
periodically extended to all other periodicity cells [13J . 
The delta function potential induces cusps in the even 



orbitals (/)^ at cc = and more generally x which are 
integral multiples of L, which are to be determined from 
4>k '(0+) - 0fc '(0-) = 2(f>k '(0+) = 2Ab^(/.fc(0) where the 
prime denotes the derivative, thus determining Ak and 
hence the normalization constant Ak- One finds Ak = 
k/\u and Nk = ^2/L{l + Ai) = Xm^2/L{\1 + k^). 
As in [12], the odd orbitals (jyjT do not see the delta 
potential and are cusp-free. 

Impurity-boson distribution function: By Eq. ([4]), the 
impurity-boson distribution function pu in the labora- 
tory system is the expectation value oiipp{x—y)ipF{x—y) 
where y is the impurity position. Then pbi{x — y) — 
X]occl'?^fe(^ — y)P where X^occ runs over the N occu- 
pied orbitals with ^ < k < {N + 1)tt/L. By Eq. (U) 
and the argument following Eq. ([5]), the wave func- 
tion in the laboratory system is <&iab(a^i, • ' ' i^M^y) — 
L^^/'^e^'^y^p{xi - y,- ■ ■ ,xn - y) where $f is the filled 
Fermi sea. All of the linear momentum resides in the 
prefactor e*''^ since the factor ^p is now translationally 
invariant. Application of the standard definition pbi{x — 
y) = N{N -1) J |$iab(2:,X2, • • • ,XN;y)\'^dx2 ' ' -dxN Ver- 
ifies the correctness of the above expression for pbi as 
a sum over density contributions from all occupied or- 
bitals. In the thermodynamic limit the sum over the 
trivial orbitals with Ak — is an integral over the Fermi 
sea from to kp, yielding §[1 — jo{2kp{x — y))] where 
jo(z) = sin(z)/2: and n = -j- = -^ is the mean fermion 
number density. In the limiting cases Xbi — Xti/kp = 
and Xbi -^ 00 the sum over the orbitals with Ak =^ 
can also be evaluated analytically in the thermodynamic 
limit; combining with the above contribution one finds 
Pbi{x — y) = kp/w = n for Xbi = as expected, and for 
Xbi -^ 00 (impenetrable impurity) one finds pbi {x — y) = 
n[l — jo(2fci?a;)] which vanishes when a; — y == as ex- 
pected. For intermediate values of Xbi the sum over the 
^^-^^ — 200 lowest orbitals with Ak ^ (!) has been eval- 
uated numerically, and adding the above expression for 
the Ak — contribution yields the results shown in Fig. 
[2]for pbi{x—y). These results differs from the first section 
(one-fermion case of the boson-impurity mixture theory 
'T'l) because the many-body wave function there repre- 
sented the impurity "dissolved" into the TG gas with 
the total linear momentum q shared equally by the im- 
purity and N bosons, whereas here the impurity moves 
through a nonmoving sea of fermions. 

Generalization to q > 0: Here we generalize to the 
case of a moving impurity by including the term — 22^ in 

Eq. ([6]). Dropping p^p/2mi and the c-number -S—, one 
finds that the A^-particle ground state is a filled Fermi 
sea of the lowest N orbitals 0^ {x) of the single-particle 
Schrodinger equation Hq(j)k{x) = {k^ /2)(j)k{x) with 
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where Q — q/mi. The momentum shift Q can be 




thermodynamic limit where N -^ oo and L — > oo with 
n — N/L fixed. The ^^ vanishes when acting on |^o)> 



FIG. 2: Impurity-boson distribution function for various 
choices of impurity-boson coupling constant Xti- Thin solid 
line: Xti = oo from analytical expression. Other curves from 
numerical calculations with ^^-^ = 200 orbitals: thick solid 
hne Xbi ~ 0.1, dashed line Xm — 1, dot-dashed line Xm = 10. 



eliminated by a gauge transformation, writing 4>k{x) = 
e''^^Mfc(a;), after which u/j is an eigenfunction of the orig- 
inal Schrodinger equation ^. However, now periodicity 
of 0fc is not automatic, but must be enforced, leading 
to nontrivial allowed values of k determined by a secular 
equation. This greatly complicates the solution. How- 
ever, simple results can be obtained in the thermody- 
namic limit by noting that the total momentum q is nec- 
essarily an integral multiple of 27r/L, i.e., q = v2-k jL with 
V a non-negative integer. For finite, nonzero g, the ther- 
modynamic limit is then to be attained by letting the in- 
teger 1/ — > oo and fermion number iV -^- cx3 in such a way 
that V jN — > q/2Tin where n = N/L is the fermion den- 
sity in the thermodynamic limit where N —f oo, L —^ oo, 
N/L -^ n. In that limit there are infinitely many values 
of rrii which are divisors of ;/ = Lq/2-K, and the set of all 
such divisors is dense on the positive real axis. For all 
such jTii, Q is an integral multiple of 2t:/ L, the orbitals 
reduce to the q = Q orbitals ([8]), the mean fermion mo- 
mentum {pf) — 0, and the impurity-boson distribution 
function reduces to that of Fig. [51 

Momentum fluctuation and thermodynamic limit: Here 

-2 

we show that the term ^^ in ^ is negligible, {pp) is 
zero since the orbitals ^ are real and have zero flux; 
hence {p^p) = (\I'o|p|-|'J'o) is the mean square fluctua- 
tion of the fermion momentum about its mean. The A'^- 
fermion ground state is |4'o) = (11^^1)10) where the c\ 
are creation operators for the N lowest orbitals 4)k of Eq. 
dH}. One has pp = Lfefe' 4(<^fclT J;l</'fe')c/c'- The diag- 
onal elements (0(+)|i^|0(+)) and (</>(-) |i^|^(-)) van- 
ish since their integrands are odd, and the off-diagonal 
elements lead to pp — X^fc 



-ikAi^ 



.g(+)t2(-) 



where 0{L^^) is a func- 



since both cf)^^' and 4>\. are occupied below the Fermi 
surface and both empty above it. There are two factors 
of L~^ from the other term and two factors of L/2tt from 
Efcfc, - iL/27:)^Jdkdk', so ^^{M>o\p%\^o) = Oil) in 
the thermodynamic limit, negligible compared with the 
total ground state energy which is 0{L). We expect that 
the contribution of this term to pbi(x — y) is only 0(L~^). 
The mean square fluctuation of the impurity momentum 
Pi is also 0(1) since the total momentum q is conserved. 

Prospects: We have obtained exact results in the 
thermodynamic limit for the ground state energy and 
impurity-boson distribution function of an impurity in 
a TG gas, for arbitrary impurity-boson mass ratio and 
interaction strength. Control of these two additional pa- 
rameters should provide an arena for future experiments 
building on the pathbreaking experiments M, |5(. 

The first part of this paper is an outgrowth of work [7[ 
initiated when we were participants in the 2007 workshop 
"Quantum Gases" at the Institut Henri Poincare-Centre 
Emile Borel (IHP) in Paris, and we are grateful to the 
workshop organizers and the IHP for hospitality and sup- 
port. We also thank Ewan Wright for helpful comments. 
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